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Abstract 

The distributions of the Yang-Lee zeros of the ferromagnetic and antiferromagnetic Q-state Potts 
models in one dimension are studied for arbitrary Q and temperature. The Yang-Lee zeros of the 
Potts antiferromagnet have been fully investigated for the first time. The distributions of the 
Yang-Lee zeros show a variety of different shapes. Some of the Yang-Lee zeros lie on the positive 
real axis even for T > 0. For the ferromagnetic model this happens only for Q < 1, while there 
exist some zeros of the antiferromagnetic model on the positive real axis both for Q < 1 and for 
Q> 1. 
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I. INTRODUCTION 

The Q-state Potts model □ DDDD is a generalization of the Ising (Q = 2) model. 
The Q-state Potts model exhibits a rich variety of critical behavior and is very fertile ground 
for the analytical and numerical investigations of first- and second-order phase transitions. 
The Potts model is also related to other outstanding problems in physics and mathematics. 
Fortuin and Kasteleyn 0, 3} have shown that the Q-state Potts model in the limit Q — > 1 
defines the problem of bond percolation. They |7[ also showed that the problem of resistor 
network is related to a Q = limit of the partition function of the Potts model. In addition, 
the zero-state Potts model describes the statistics of treelike percolation & and is equivalent 
to the undirected Abelian sandpile model M. The Q = \ state Potts model has a connection 

n 

to a dilute spin glass [lflj. The Q-state Potts model with < Q < 1 describes transitions in 



The partition function 
or the Whitney 



le pe 

Q 



the gelation and vulcanization processes of branched polymers 
of the Potts model is also known as the Tutte dichromatic polynomial 
rank function in graph theory and combinatorics of mathematics. 

By introducing the concept of the zeros of the partition function in the complex magnetic- 
field plane, Yang and Lee [14J proposed a mechanism for the occurrence of phase transitions 
in the thermodynamic limit and yielded a new insight into the unsolved problem of the 



ferromagnetic Ising model in an ar 
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jitr ary nonzero external magnetic field. It has been 
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23, 



24j that the distribution of the zeros of a model 
determines its critical behavior. Lee and Yang also formulated the celebrated circle 
theorem which states that the partition function zeros of the Ising ferromagnet lie on the 
unit circle in the complex magnetic-field (x = e^ H ) plane. However, for the ferromagnetic 
Q-state Potts model with Q > 2 the Yang-Lee zeros in the complex x plane lie close to, 
but not on, the unit circle with the two exceptions of the critical point x = 1 (H = 0) 



itself and the zeros in the limit T = 
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It has been shown 



27, m 



that the distributions of the ferromagnetic Yang-Lee zeros for Q > 1 have similar properties 
independent of dimension. 

Recently, the exact results on the Yang-Lee zeros of the ferromagnetic Potts model have 
been found using the one- dimensional model. Mittag and Stephen |29| studied the Yang-Lee 
zeros of the three-state Potts ferromagnet in one dimension. Glumac and Uzelac Js(| found 
the eigenvalues of the transfer matrix of the one-dimensional Potts model for general Q. 
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In particular, they have shown that the ferromagnetic Yang-Lee zeros can lie on the real 
axis for Q < 1. However, the properties of the Yang-Lee zeros of the antiferromagnetic 
3 otts model have never been known except for the one- dimensional Ising antiferromagnet 
31( . In this paper, we study the ferromagnetic and antiferromagnetic Yang-Lee zeros of the 
one-dimensional Q-state Potts model. 

II. PARTITION FUNCTION 

The Q-state Potts model in an external magnetic field H q on a lattice G with N s sites 
and Nf, bonds is defined by the Hamiltonian 

= -Jj2 S ( a i> a i) ~ H q J2 S ( a k,q), (1) 

(id) k 

where J is the coupling constant (ferromagnetic model for J > and antiferromagnetic 
model for J < 0), indicates a sum over nearest-neighbor pairs, ai = 1, 2, Q, and q is 
a fixed integer between 1 and Q. The partition function of the model is 

Zq=Y. ^ (2) 

where {cr n } denotes a sum over Q Ns possible spin configurations and (3 = (fc B T) _1 . The 
partition function can be written as 

Z{a, x,Q) = Y,Y, n Q( E ' M)a E x M , (3) 

E=0 M=0 

where a = y" 1 = e^ J , x = e /3Hq , E and M are positive integers < E < iV& and < M < N s , 
respectively, and Qq(E,M) is the number of states with fixed E and fixed M. The states 
with E = (E = Nf,) correspond to the antiferromagnetic (ferromagnetic) ground states. For 
ferromagnetic interaction J > the physical interval is 1 < a < oo (oo > T > 0), whereas 
for antiferromagnetic interaction J < the physical interval 0<a<l(0<T< oo). The 
parameter Q enters the Potts model as an integer. However, the study of the Q-state Potts 
model has been extended to continuous Q due to the Fortuin-Kasteleyn representation of 
the partition function (| ?1 and its extension j^ . 

For the one- dimensional Potts model in an external field the eigenvalues of the transfer 
matrix were found by Glumac and Uzelac j^. The eigenvalues are \± = (A±iB)/2, where 
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A = a(l + x)+Q-2 and B = -iy/[a(l - x) + Q - 2] 2 + 4(Q - l)x, and A = a - 1 which is 
(Q— 2)-fold degenerate. The partition function of the one-dimensional model (N = N s = N b ) 
is given by 

Z N (a,x,Q) = X^ + X N + (Q-2)X^. (4) 



III. PARTITION FUNCTION ZEROS 



When A + and A_ are two dominant eigenvalues, the partition function becomes 



(5) 



for large N. If we define A = 2C cosip and B = 2C smip, where C = (a — l)(a + Q — l)x, 
then X± = C exp(±iip), and the partition function is 



Z N = 2C N cos Ni/j. 

The zeros of the partition function are then given by 

2k + 1 



(6) 



2N 



-7T, 



k = 0,l,2,...,iV- I. 



(7) 



In the thermodynamic limit the locus of the partition function zeros is determined by the 
solution of 

A = 2Ccosij, (8) 

where < ip < ir. In the special case Q = 2 the contribution by the eigenvalue Ao disappears 
from the partition function, Eq. (4), and the equation (8) determines all the locus even for 
finite systems. From Eq. (8) the locus of the Yang-Lee zeros for any Q is obtained to be 



xi(il>) = 



1 



ficosijj + iJ f 2 



(9) 



where fi = (a — l)(a + Q — 1) and f 2 = fi sin 2 tp + Q — l. The edge zeros of xi(ip) are given 
by 

1 



/(o-i)( fl + Q-i)±^g 

from Xi(0) and Xi(n). If fi > and f 2 > or f\ < and / 2 < 0, it is easily verified that 

a + Q - 2 



(10) 



(ii) 
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From Eq. (11) we see that the zeros of x\(ip) lie on a circle in the complex x plane. The one 
point of the circle xi(ip), 

/tt\ a + Q-2 

Xl UJ = — — (12) 

always lies on the real axis. However, if f\ > and f% < or A < and f'2 > 0, the zeros 
of Xi(ijj) lie on the real axis. Recently, Ghulghazaryan et al. |33J] tried to understand the 
properties of Eq. (9) in the antiferromganetic (a < 1) region. 

On the other hand, when A + and Ao are two dominant eigenvalues, the partition function 
can be written as 

Z N ~ A^ + (Q - 2)A^ (13) 
for large N. The partition function zeros are then determined by 

A± = (2-Q)Wex pW >), (14) 
^0 

where 

= 0fe = 7T' k = 0,l,2,...,N-l. (15) 
In the thermodynamic limit the locus of the partition function zeros is determined by the 
solution of 

aV + {Q- l)x - axA + (a - l)Ae i4> - (a - l) 2 e 2i(f> = 0, (16) 

where < < 2tt. The equation (16) also determines the locus of the zeros when A_ and 
Ao are two dominant eigenvalues. Eq. (16) gives the second locus of the Yang-Lee zeros 

e*+[a + Q - 2 - (a - l)e i(t> ] 
a + Q — 1 — ae 1 

which is the unit circle in the limit a — > 00 and a type of limagon of Pascal for a = 0. The 
two points of the locus X2 (</>), 

z a (0) = 1 (18) 

and 

2a + g-3 

X2(7F) = ~ 2a + Q-V (19) 
can be lie on the real axis. From |A±| = |Ao| we also obtain 

a — 1 . 

x, = — -. 20 

a+Q-1 ( J 

The second locus, Eq. (17), of the Yang-Lee zeros, although it is very important in under- 
standing the Q-state Potts model, has never been considered in the literature until now. 



= 1 ' ( 17 ) 



IV. FERROMAGNETIC (a > 1) YANG-LEE ZEROS 

In this section we study the Yang-Lee zeros of the ferromagnetic Q-state Potts model for 
a > 1. 



A. Q > 1 

In this case fi and /2 are always positive and the Yang-Lee zeros Xi(i/j) lie on a circle 
where the eigenvalues satisfy 



I A. 



|Ao| \ 



(a + Q-l)(a + Q-2) 
a(a — 1) 



> 1, 



(21) 



which implies that the locus £2(0) does not appear. For 1 < Q < 2 the zeros lie inside the 



unit circle while for Q > 2 the zeros lie outside the unit circle, as shown in 
special case Q = 2 we of course find that = 1, as proved by Lee and Yang 

The argument 9 of xi(ip) (= \xi(ip)\e ie ^) is given by 



23, |30(Lln the 

la. 



e 

cos - 



(a - l)(a + Q-l) 



cos i/j. 



(22) 



2 ^ a(a + Q - 2) 

As a — > 00 (T — > 0), the zeros approach the unit circle. For a = 1 (T = 00), cos(#/2) = 0, 
so #(V0 = 7T and all the zeros lie at 1 — Q. As Q — > 00, the radius |xi( , ?/')| increases without 
bound. The Yang-Lee edge zero for T > is given by 



#n = 2 cos" 1 



(o-l)(a + Q-l) 



> 0, 



(23) 



\ a(a + Q-2) 

while 6q = at T = 0. Therefore, we conclude that no zero lie on the positive real axis for 
any T > 0. 



B. Q < 1 



At a = 1 (T = 00), all the Yang-Lee zeros lie at the point 1 — Q on the positive real 



axis between X- and x + 



axis. For l<a<l — ^ + all the zeros lie on the positive rea 
(0 < x_ < x + < 1), as pointed out by Glumac and Uzelac [30L 133. 
the locus consists of the loop X2{4>) (0* < < 27r — 0*) and the line x\{jp) (0 < ip < -0*) on 



3j. For a > 1-f + 
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the positive real axis between x* = xi(-0*) and x + = xi(0) (0 < x* < x + < 1). The loop x 2 
meets with the line x\ at the point 

x* = x 2 (0*) = x 2 (2ir - 0*), (24) 

where |A+| = |A_| = | Ao j - The loop x 2 (0) cuts the real axis at two points x 2 {tx) and x* 
(> x 2 (n)). The sign of x 2 (rr) is positive forl-^ + ^<a< ^ and negative for a > 
Figure 1 shows the locus of the Yang-Lee zeros for Q = \ and a = | from which we obtain 
x 2 (ti) = yg, x* = |, and x + = 0.8119. 

At Q = 0, the line Xi(ip) shrink to the point 

x* = x + = x 2 (0) = 1, (25) 

and all the zeros lie on the loop x 2 (0) (0 < < 2rr) for a > 1. As Q — > 1, two points x* and 
x + approach (a — l)/a together, and the line xi(-0) disappears. At a = oo (T = 0), the line 
Xi(ip) again shrink to the point x* = x + = 1, and the loop x 2 (0) becomes the unit circle 
centered at the origin. 

V. ANTIFERROMAGNETIC (0 < a < 1) YANG-LEE ZEROS 

In this section we investigate the Yang-Lee zeros of the antiferromagnetic Q-state Potts 
model for < a < 1. 

A. Q > 1 

Because fi is always negative, the zeros of Xi(ip) lie on a circle if f 2 < (-00 < "0 < t 1 " - "0o) 
and on the real axis if f 2 > (0 < -0 < ip or ir — ip < -0 < it). For a < 1 — ^, the locus 
consists of the line Xi(-0) (0 < -0 < -0 O and 7r — -0 O < -0 < ip*) on the negative real axis 
between x + = Xi(0) and x* = Xi(-0*) (x + < x* < 0), the circle Xi(-0) (V'o < "0 < tt — -0o) 
with the radius |xi(-0)|, and the loop x 2 (0) (0 < < 0* and 2n — 0* < < 27r), inside the 
circle xi. The loop x 2 again meets with the line x\ at the point 

X* = X 2 (0*) = X 2 (27T - 0*), (26) 
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where |A+| = |A_| = |A |. The circle cuts the real axis at two points x? and xi(|) (= — x\ > 
1), where the point is defined by 

x° = xi(^b) = x 1 (ir- Vo) = a + Q ~ 2 . (27) 

Similarly, the loop also cuts the real axis at two points x* (< 0) and x 2 (0) (=1). At a = 
(T = 0), the circle xiijp) disappears for 1 < Q < 2, and the locus consists of the line x\(ip) 
on the real axis between — oo and x* (< 0) and the loop x 2 (0). For example, figure 2 shows 
the locus for Q = ^ and a = j^. In this case we obtain x + = —54.83, x\ = —8, x* = — |, 
the radius \x x \ = 8 for the circle, = 48.19°, ^* = 134.07°, and 0* = 46.02°. 
At a — 1 — ^, two points x* and x° on the real axis meet at 

= x? = -1, (28) 

other two points Xi(|) and £2(0) on the real axis also meet at 

*iQ=s 2 (0) = l, ( 29 ) 

and the circle x\{jp) and the loop #2(0) become the identical locus as the unit circle. On 
this unit circle, three eigenvalues have the same magnitude 

|A+| = |A_| = |A | = |. (30) 

Therefore, the locus consists of the line x\(ip) on the real axis between x + and x* and the 
circle. 

In the region 1 — Q- < a < I — ^ + the circle Xi(-0) disappears, and the line Xi(ip) 
(0 < ip < ip*) on the negative real axis between x + and x* (x+ < x* < 0) again meets with 
the loop #2(0) (0* < < 27r — 0*) at the point x*. The loop cuts the real axis at two points 
x* and 2:2(71"). The sign of x 2 (tt) is positive (0 < x 2 {n) < 1) for 1 — ^ < a < and negative 
for < a <i_<2-|_-y^2. Figure 3 shows the locus of the Yang-Lee zeros for Q = 3 and 
a = jq which give x 2 (7r) = — -y, x* = — |, and x + = —777.84. At a = (T = 0), the locus 
still consists of the line Xi(ip) on the real axis between —00 and x* (< 0) and the loop x 2 (0) 
for 2 < Q < 4. At T = 0, the sign of x 2 (7r) is positive (0 < x 2 (tt) < 1) for 2 < Q < 3 and 
negative for 3 < Q < 4. As Q — > 1, two points x + and x* approach (a — l)/a together, and 
the line xi (■?/>) disappears. 

At a = 1 — ^ + ^, x* = x_, and the loop x 2 (0) shrinks to the point x*. For a > 1 — 
the loop disappears, and the only locus is the line Xi(ip) on the negative real axis between 
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x + and x_ (x + < x_ < 0). As Q — > oo, x± —> — oo for a < 1. At a = 1 (T = oo), two edge 
zeros x+ and x_ meet, and the line xx(ifi) shrinks to the point 

x + = X- = 1 — Q (31) 

for Q > 1. In the special case Q — 2, fi < and /2 > 0, and the line xi(^) is the only 
locus of the Yang-Lee zeros. Therefore, all the Yang-Lee zeros of the antiferromagnetic Ising 
model lie on the negative real axis between x + and X- (> x + ) for a < 1, as shown by Yang 
31(. At a = (T = 0), all the zeros lie on the negative real axis between two edge zeros, 

x a = -oo (32) 

and 

Xb = {Xb ~ 0) ' (33) 

for Q = 2 and Q > 4. 



B. Q < 1 

For 0<a<l — ^ — the locus consists of the loop 22(0) (0 < (fi < 0* and 27r — 0* < 
(fi < 2n) and the line xi(ip) (0 < ^ < ip*) on the positive real axis between = Xx(ip*) and 
x + = xi(0) (1 < < x + ). The loop X2 cuts the real axis at two points £2(0) (= 1) and X* 
(> 1). The loop again meets with the line x\ at the point 

= x 2 ((fi*) = x 2 (2tc - 0*), (34) 

where |A + | = |A_| = |A |. At Q = 0, two points x 2 (0) and x* meet, the loop x 2 ((fi) shrinks 
to the point 22(0) = and all the Yang-Lee zeros lie on the positive real axis between x_ 
(= x* = 1) and x + (> X-) for a < 1. 

At a = 1 — y — = x + , and the line Xi(ip) shrinks to the point x*. In the 

region 1 — y — ^ < a < 1 — the line disappears, and the locations of Yang-Lee zeros 
are completely determined by the locus X2(0) (0 < (fi < 2tt) which cuts the real axis at 
two points X2(0) (= 1) and X2(7r) (> 1). The shape of the locus X2{4>) is changed from a 
waterdrop-like shape for 1 — ^ — ^ < a < 1 — Q through a circle with center arid 
radius ^zq for a = 1 — Q to a crescent-like shape for 1 — Q < a < 1 — 2 ■ Fig ure 4 shows 
the locus of the Yang-Lee zeros for Q = and a = | which give X2(7r) = y. 
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At a — 1 — ^, the loci Xi(ip) and x 2 (0) become identical to be the unit circle on which 
three eigenvalues again have the same magnitude 

|A+| = |A_| = |A | = |. (35) 

For a > 1 — ^, all the Yang-Lee zeros lie on the circle Xi(ip) (0 < ip < it) which cuts the 
positive real axis at the point a?i(f) (0 < Xi(§) < 1)- At a = 1 (T = oo), this circle shrinks 
to the point Xi(ip) = 1 — Q for Q < 1. 

VI. ONE-STATE (Q=l) POTTS MODEL 

At Q = 1, the partition function becomes 

Z N (a,x,Q = 1) = (ax) N . (36) 

Therefore, all the Yang-Lee zeros lie on the point x — 0. On the other hand, in the limit 
Q — > 1, the line xi(^) disappears, and the circle x\{ip) is the uniform distribution of the 
zeros because the circle approaches 

a — 1 



xi(ip) 

The loop x 2 (0) also approaches 



e 2 T (37) 



z 2 (0) = — e*, (38) 



which is the same as the circle xi (■?/>). 



VII. CONCLUSION 



We have studied the interesting properties of the ferromagnetic and antiferromagnetic 
Yang-Lee zeros of the one-dimensional Q-state Potts model for arbitrary Q and temperature. 
It has been shown that the distributions of the Yang-Lee zeros have a variety of different 
shapes. In particular, the antiferromagnetic Yang-Lee zeros of the Potts model have been 
fully investigated for the first time. One of the most interesting results is that some of 
the Yang-Lee zeros lie on the positive real axis even for T > 0. For the ferromagnetic 
Yang-Lee zeros this happens only for Q < 1, while there exist some Yang-Lee zeros of the 
antiferromagnetic Potts model on the positive real axis both for Q < 1 and for Q > 1. 
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The results obtained from the one-dimensional Potts model may be considered as a road 
to the full understanding of the Yang-Lee zeros in higher dimensions in that for Q > 1 the 
distributions of the ferromagnetic Yang-Lee zeros have similar properties independent of 
dimension j^i 28^ . 
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FIG. 1: The locus of the Yang-Lee zeros in the complex x = e@ Hq plane for Q = ^ and a = e@ J = |. 
For comparison, the zeros for a finite-size system (N = 100) are also shown (open circles). 
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FIG. 2: The locus of the Yang-Lee zeros for Q = and a = j^. 
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